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■ Abstract 

> 

' We analyze the structure of the vacuum polarization tensor in the presence of a background 

. electromagnetic field in a medium. We use various discrete symmetries and crossing symmetry to 

1/^ ' constrain the form factors obtained for the most general case. From these symmetry arguments, we 

, show why the vacuum polarization tensor has to be even in the background field when there is no 

On ' background medium. Taking then the background field to be purely magnetic, we evaluate the vacuum 
polarization to linear order in it. The result shows the phenomenon of Faraday rotation, i.e., the 

' rotation of the plane of polarization of a plane polarized light passing through this background. We 

l| find that the usual expression for Faraday rotation, which is derived for a non-degenerate plasma in 

^ ^1 the non-relativistic approximation, undergoes substantial modification if the background is degenerate 

^ ' and/or relativistic. We give explicit expressions for Faraday rotation in completely degenerate and 

, ultra-relativistic media. 



^ ' 1 Scope and outline of the paper 
H ; 

It is well known that electromagnetic wave propagating through a medium in an ambient magnetic 
field suffers Faraday rotation, i.e., the plane of a plane polarized light rotates as it travels through the 
medium in the magnetic field. The amount of this rotation is derived in various texts on electromagnetic 
theory |Q] and plasma physics Q with the assumptions that the medium consists of non-relativistic and 
non-degenerate electrons and nucleons. 

Faraday rotation is extensively used in a variety of situations, including astrophysical and cosmological 
ones 1^, |3|. In such situations, either of the aforesaid assumptions about the medium may not be valid. 
For example, for compact stars, the plasma is likely to be degenerate. In the very early universe, when 
the temperature was high, the assumption of non-relativistic plasma is bound to break down. Motivated 
by such situations, we reinvestigate this problem. For a general framework, the formalism of quantum 
field theory proves to be helpful. The aim of this paper is to use quantum field theoretical formalism to 
calculate Faraday rotation in different kinds of media. 

The paper is organized as follows. In Sec. ^ we introduce the vacuum polarization tensor and find 
its most general form in a background medium and in the presence of a general electromagnetic field, 
consistent with Lorentz and gauge invariances. Calculation of the vacuum polarization tensor requires 
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the electron propagator, which is discussed in Sec. ^. We summarize there how Schwinger's proper-time 
propagator in a constant magnetic field is modified in the presence of a background medium. In Sec. ^, 
we show how some discrete symmetries help us constraining some of the form factors appearing in the 
vacuum polarization. Following this, we set up the calculation in Sec. |^. Starting from the basic Feynman 
rules, we arrive at an expression for the polarization tensor that is explicitly gauge invariant. In Sec. ^ we 
derive the expression for Faraday rotation per unit length in terms of the components of the polarization 
tensor. Then in Sec. |^, we provide explicit results for different kinds of backgrounds. This is the section 
which contains the essential results of the paper. The non-relativistic and non-degenerate case is shown 
in Sec. |7.3| , where we obtain the result usually quoted in textbooks. In Sees. 7_A and 7^, we find results 
for a completely degenerate medium and an ultra-relativistic one. Finally, we present our conclusions. 



2 Form factors in the polarization tensor 

The classical action of a free electromagnetic field is given by 



where 



A^-^Jd^xF^p{x)F^P{x). (2.1) 
In the momentum space, this can be written as 

where C is the momentum-space Lagrangian, which can be obtained by taking Fourier transforms in 
Eq. (13): 

C = ~h'gx,A\k)AP{-k), (2.3) 

5Ap = 5Ap-^. (2.4) 

Once quantum corrections are added, one obtains more quadratic terms in the Lagrangian. These are 
represented by the vacuum polarization tensor Il\p. In other words, after the quantum corrections are 
put in, the quadratic part of the Lagrangian becomes 

£ = i [-k^gxp + U^p{k)] A\k)AP{~k) . (2.5) 

Owing to gauge invariance, II^p satisfies the conditions 

fc^nAp(fc)-o, KnAp(fc) = o. (2.6) 

In addition, Bosc symmetry implies 

nAp(fc) = npA(-fc). (2.7) 

In the vacuum, the tensor II^p depends only on the momentum vector k. Thus, the most general form 
for Hxp is given by 

Hap = Ho [k^gxp - kxkp] , (2.8) 

where IIq is a Lorentz-invariant form factor, which can therefore depend only on fc^. The important 
point is that the tensor structure for IIap is exactly the same as the tensor appearing in the classical 
Lagrangian. Thus, this correction term can be done away with by a redefinition of the photon field A'*'. 
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In a nontrivial background, this is no more the case. Although Il\p still has to satisfy Eq. (|2£ 
the form given in Eq. ( ^.8D does not follow. This is because Hxp can now depend, apart from the 
momentum vector k^, on various vectors or tensors which characterize the background medium. Even 
for a homogeneous and isotropic medium, there is an extra vector in the form of the velocity of its center 
of mass, u^. The most general form for IIap in the presence of these two vectors has been discussed in 
the literature [Q. 

Our interest lies in a more complicated background where in addition to the medium there is also an 
external electromagnetic field Bxp. We will work in the weak field limit throughout. This means that the 
background field will be considered feeble, and we will keep only linear terms in it. For the moment, we 
will not specialize to magnetic fields. We will keep the discussion general, with a medium characterized 
by the vector and a background electromagnetic field B\p in arbitrary direction. 

There are many independent tensors constructed out of fc'*', and Bxp which satisfy Eq. (2^). For 
future convenience, we categorize them into several groups. In the first group, there is only one tensor 
which depends only on the vector k'^, viz., 

Pxp =9Xp = 9Xp-^- (2.9) 



This, of course, is the same tensor which appears in Eq. (|2.8| ). 

In the second group, we include tensors constructed only of k and u. These are Q: 

PxJ = «AS^p/^^^ (2.10) 

P[^J = expark^u-", (2.11) 

where 

ux=gxau''. (2.12) 

Next, we bring in tensors constructed from k and B only, without any occurence of u. Using the 
shorthand 



we can write these as: 



ik-B)x - fc^B.A, (2.13) 



- k^Bxp^kx{k-B)p + kp{k-B)x, (2.14) 



xp ^ -°Ap - i^xyi^ ■ o)p 

y(2 
\p 

One might think that there might be additional terms obtained by replacing Bxp by -Bap, where 



P'Z^ = expark^k-BY . (2.15) 



-Bap — -SXparB'^'^ . (2-16) 

But it is straight forward to show that no other independent term arises this way. 

Finally, to write down the tensors where all three of fc, u and B occur, we employ a notation (u ■ B)x 
defined in a way similar to (fc • B)x- Then the tensors are: 

P^'j') = k-uBxp-ux[k-B)p + Up{k-B)x (2.17) 

P'lf^ = expark''{u-By (2.18) 

= ux{k-B)p-Up{k-B)x (2.19) 

Pxf^ = ux{k-B)p + Up{k-B)x (2.20) 

P'lf = ux9pr{u-BY -Up-gxAu-Bf (2.21) 



3 



p"(6) 


= uxgpriu-B] 


^ + Upfjxriu 


BY 


p"(7) 
^Ap 


= uxik-B)p- 


- Upik - B)x 




p"(8) 
^Ap 


= UA(fc-S)p4 


- Upik - B)x 




p"(9) 
^Ap 


= Uxgpriu-B] 


^ - Upgxriu 


BY 


p//(10) 


= Uxgpriu-B] 


^ + Upfjxriu 


BY 



(2.22) 
(2.23) 
(2.24) 
(2.25) 
(2.26) 

A collection of gauge invariant tensors which can appear in the vacuum polarization were listed by 
Perez Rojas and Shabad ||^. They have tensors involving more than one powers of the background field, 
which we are not interested in. As for the tensors linear in Bxp, they list what we call Pxl\ P')}iP 



P'xp \ but none of the rest. We conclude that the most general form for IIap consistent with Eq. (2.6) 



and keeping only terms linear in the external electromagnetic field is given by: 

n,, = n(')p« + Y: ^'""p^ + E n"^M?^ , (2.27) 

i i i 

where in each case, the sum over i runs over the appropriate set of values. The coefficients of the tensors 
are form factors, which we discuss now. 

First, notice that the tensors P^J and P^^^J do not depend on the background electromagnetic field. 
The form factors associated with these terms are related to the dielectric constant and the magnetic 
permeability of the medium |^ . The tensor pj^J also does not involve the background electromagnetic 
field. It is, however, parity asymmetric, and accounts for natural optical activity [Q. Since our aim is 
to discuss Faraday rotation which is also a type of optical activity, we will disregard any natural optical 
activity. Thus, we assume that the form factor associated with P^J is zero for our medium. 

The form factors are Lorentz invariant quantities. Thus, they can depend only on the Lorentz invariant 
combinations of fc^, u^ and Bxp. Since u^ — 1, we can obtain the following invariant parameters, keeping 
at most one factor of the background field: 

UJ = k-u, (2.28) 
K = ^JJ^ - fc2 , (2.29) 
h = k^uPBxp, (2.30) 
b = k^uPBxp- (2.31) 

In addition, of course, the form factors can depend on the Lorentz scalars which define the background 
medium, e.g., the chemical potential ^ and the temperature 1/(3. 

Since we are interested only in linear terms in the background field, we have discarded higher order 
invariants in Bp^. Moreover, notice that the tensors P'^^^ and P'xp^^ are linear in the background field. 
Thus, for their co-efficients, we can neglect the field dependence for the sake of consistency. Thus, for our 
purpose, the form factors 11'^*^ and n"*^*) should be treated as functions of uj and K only, and possibly of 
H and p. We summarize this statement as: 

n« =n«(c.,if,^,/3), (2.32) 

where Yl'-'^ stands for the n'(^)'s and the n"(^)'s. In Sec. 0, we will see how arguments about various 
discrete symmetries restrict the form factors in significant ways. The constraint of hermiticity of the 
action implies 

nAp(fc)-n;,(fc), (2.33) 

whose consequences will be mentioned at the end of Sec |[ 
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Figure 1: One- loop diagram for vacuum polarization. 

3 The electron propagator 

At the 1-loop level, the vacuum polarization tensor arises from the diagram in Fig. |l|. The dominant 
contribution to the vacuum polarization comes from the electron line in the loop. To evaluate this diagram, 
one needs to use the electron propagator within a thermal medium in the presence of a background 
electromagnetic field. Rather than working with the complicated expression for a general background 
field, we will specialize to the case of a purely magnetic field. Once this is assumed, the field can be taken 
in the z-direction without any further loss of generality. We will denote the magnitude of this field by B. 

Ignoring at first the presence of the medium, the electron propagator in such a field can be written 
down following Schwinger's approach [0, ||, 

iSlip)^ / rfse*(^''^)Gb,s), (3.1) 
Jo 

where $ and G are defined below. To write these in a compact notation, we decompose the metric tensor 
into two parts: 

ff"/3 = gin - gap , (3-2) 

where 

glf, = diag(l, 0,0,-1), 

gif^ = diag(0, 1,1,0). (3.3) 
This allows us to write, for any two objects q and q' (including the 7-matrices) carrying Lorentz indices, 

Q-q'w = %% - Qsl's : (3-4) 

q-q'± = QiQ'i + 9292 ■ (3-5) 

Using these notations, we can write 

Hp,s) ^ ^s(pf^-^-^^^pl-m'^-e\s\, (3.6) 

^ieSsaz / ^ — ieBsGz 



G{p,.s) = ——^ U 



where 



and we have used 



cos(eSs) \ " cos(eSs) 
= {I + i(Tzta.neBs){-^\\ + m) — {sec^ eBs)^± , (3.7) 



«7i72 = -707375 , (3.8) 



jeBsa, = COS cBs + ia, siu eBs . (3.9) 
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The expression for $ can have an additional gauge dependent phase factor, but it does not contribute 
to the polarization tensor. Usually, one writes s instead of \s\ in Eq. ( |3.6| ). It is equivalent since in the 
range of integration indicated in Eq. s is never negative. However, the definition of s) is useful 
in this form for what follows next. 

In the presence of a background medium, the above propagator is modified to pO|: 



iS{p) ^ iSsip) ~ riF{p) iSB{p)~iSB (p) 



(3.10) 



where 



5B(p)^7o4b)7o (3.11) 

for a fermion propagator, and "qpip) contains the distribution function for particles and antiparticles: 

t^f{p) = e(p • u)fF{p. M,/3) + Q{-P ■ u)fFhP, -/i,/3) . (3.12) 

Here, 6 is the step function, which takes the value +1 for positive values of its argument and vanishes 
for negative values of the argument, and fp denotes the Fermi-Dirac distribution function: 

1 



fF{p,^^,f3) 



(3.13) 



Putting in the form of Sg(j)) from Eq. (3.1), we obtain the additional term in the propagator to be 



Slip) ^ -^VF{p)[S^{p)~sl{p) 

/oo 
ds e*(f'^)G(p,s), 
-OO 



(3.14) 



with $(p,s) and G{p,s) defined in Eqs. (M) and (3.7) 



It is straight forward to see that when B = 0, the propagator in Eq. (3.1) reduces to 

/•CSC 

iSq {p) = I ds exp [is (p^ — + ie)] (^ + rn) 
Jo 



^ + TO 



p^ — m? + ie 



which is the vacuum propagator. In the same limit, the background dependent part reduces to 

S';'(p) = -27r(5(p2 _ m')i^F{p){i> + ™) • 



(3.15) 



(3.16) 



4 Discrete symmetries and the form factors 

Before embarking on diagram calculations, let us discuss some of the symmetries of the problem, which 
will help us constraining some of the form factors. 



4.1 Bose symmetry 

This symmetry was already discussed in Eq. ( ^.7] ), viz., under the operation 

k — fc , A <-+ p . 

the vacuum polarization tensor must be invariant. 
y(0 



(4.1) 



Since the tensors P^^ change sign under the operation of Eq. (4.1), this impHes that the associated 
form factors should satisfy the condition 



H'W (c., if , /i, /3) = - H'« (~a;, if , /3) . 
6 



(4.2) 



On the form factors denoted by 11", the effect is more comphcated since some of the tensors P"^ are 
symmetric under the operation of Eq. (2.7) and some are antisymmetric. In general, let us write 



where, by inspection, we see that 



+1 for i = 1,2,3,6,7,10 

-1 for i 4,5,8,9. 



(4.3) 



(4.4) 



The associated form factors should then satisfy the relation 

n"(^) {lj, K, fi, 13) = n, n"(') {~Lo, K, /i, P) . (4.5) 

To apply this symmetry on the unprimed form factors, one has to take into account the dependence 
of these form factors on b and h defined in Eq. ( 2.31 ). Thus, this symmetry implies 

n(')(w,i^,Ai,^,6,6) = +n(')(-w,i^,Ai,/3,-6,-6) fori = 0,1 (4.6) 
n(2)(c^,X,^,/?,6,6) = -n(2)(-c^,if,^,/3,-6,-6). (4.7) 



4.2 Charge conjugation symmetry 

In calculating the form factors, we are neglecting any corrections coming from weak interactions. In 
fact, these corrections occur only at the 2-loop level, and therefore are anyway irrelevant for the 1-loop 
calculation that we will be performing. In this case, the interactions are all purely electromagnetic, and 
so they obey charge conjugation (or C) symmetry. The conclusion of this symmetry is that, lixp should 
be invariant under the substitutions: 

H^-H, B^r^-Bar- (4.8) 

Said in words, it means that if we calculate the vacuum polarization in a medium with a certain back- 
ground field, it should be the same as that obtained in a charged conjugated medium with an opposite 
background field. This means that, in the vacuum polarization, there are terms even in the background 
field and even in ^, or odd in both. The terms linear in B, which we will calculate, should therefore be 
odd in /i. This implies that the primed form factors, which are independent of the background medium, 
i.e., do not contain fi, must vanish. This is known from the direct calculations of the polarization tensor 
in absence of a medium ^ pT[ . 



4.3 A symmetry of the propagator 

Lastly, notice that in the calculation of Tl\p, the center-of-mass velocity and the chemical potential 
fjL can enter only through the function rjp in the propagator. Further, from Eq. (3.12), notice that rjp is 
invariant under the following transformation: 

M, /i^— /i. (4-9) 

So the vacuum polarization must also obey this symmetry. 

For the form factors, this fact has an interesting consequence. Some of the tensors -P^'j*'* are even in 
u, some are odd. Accordingly, the form factors would satisfy 

n"(^) (a;, K, ^^, (3) = < n"(') {-lo, K, /?) , (4.10) 

where 

-1 forz= 1,2,3,4,7,8 

fl fori = 5, 6, 9, 10. ^ ■ ^ 
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Using the consequence of C-symmetry, we can conclude that, for terms hnear in B, 



n"(') (c., K, fi, P) = -<n"(') i-Lo, K, ^1, P) . (4.12) 
In order that this is consistent with Eq. ( |4.5| ), we must have 

(no sum on i). (4.13) 



Using Eqs. ( |4^ ) and ( f4.1l| ), then, we obtain that the doubly-primed form factors U"^^^ and 

n"(io) vanish. 

The analysis performed in this section is valid for a general electromagnetic fields in the weak limit. 
Although we have used the propagator in the presence of a purely magnetic field, substituting the more 
general form does not affect the arguments. We can now discuss how this analysis can be simplified 
if the background field is a purely magnetic field in the rest frame of the background medium, i.e., if 
u'^Bcrr = 0. Among the field-dependent tensors, Px'f \ P')}p'^ and -P^p^"* vanish in this case. Therefore, 
in the final count, we will have only four form factors associated with the field-dependent tensors, viz., 
n"(i) , n"*^^^ , n"'^^^ and n"'^^^ . since all of these are antisymmetric tensors, the hermiticity condition of 



Eq. (2.33) implies that the corresponding form factors must be purely imaginary in the dispersive part. 

In addition, of course, we can have the form factors 11*^°) and 11^^^ whereas II'^' vanishes because of 
our assumption of vanishing natural optical activity. 

5 Calculation of the 1-loop vacuum polarization 

5.1 Identifying the relevant terms 

The amplitude of the 1-loop diagram of Fig.|l| can be written as 

in^,{k) =-J ^(*e)2 tr [j>,tSiphpiS{p + k)] , (5.1) 
where the minus sign on the right side is for a closed fermion loop, and S{p) is the propagator given in 



Eq. (P.IOD . This imphes 



nA,(fc) = -le^ I ^ tr [7A z^(p)7p iS{p + k)] . (5.2) 



From Eq. ( |3.10 ), we see that there are two terms in the propagator — the vacuum part Sg{p) and 



the other part which involves the background matter distribution. If we insert two such propagators in 



Eq. (5.2), we will obtain four terms. 

The term obtained from the factor in both propagators is the contribution in the vacuum. It has 
no importance to our discussion of background effects. The terms with the distribution function factor 
from both propagators contributes only to the absorptive part of the vacuum polarization, which we do 
not discuss in this article. Thus we are left with the terms in which we use the vacuum part of one 
propagator and the background dependent part of the other. These terms contribute to the part 11"^ in 



the notation of Sec. ||. Thus 



n';(,(fc) = -^e' y tr [J^Sl{ph,^S^{p') + 7a ^S^iph, Sl{p')\ , (5.3) 

where, for the sake of notational simplicity, we have used 

p'=p + k. (5.4) 
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Substituting p by — p' in the second term and using the cyclic property of traces, we can write Eq. (5.3) 
as 



2 I d^V 



tr 



(5.5) 



Using now the form of the propagators from Eqs. (3.1) and (3.14), we obtain 



n^p(fc) 



d^p 



ip) tr(7AG(p,s)7pG(p',s')) + riF{-p) tr(^7pG(-p, s)7aG(-p', s') 



(5.6) 



5.2 Extracting the gauge invariant piece 

In order to discuss Faraday effect, we need only the terms in the vacuum polarization tensor which are 



odd in B. Notice that the phase factors appearing in Eq. (5.£) are even in B. Thus, we need only the 
odd terms from the traces. Performing the traces is straight forward, and the odd terms come out to be 



(M 



[p) / dse*(P'^) / ds' e*(f''^')i?A, 



where we have introduced the notation 



and 



Rxp — 



V-{p) = Vf{p) " Vf{-p) , 

£Ap03™^ ( tan eBs — tan eBs') 
+£Apa||/3|| taneSs - p'""/" taneBs' 

+e Apa II /3i (p" " P''^^ tan eSs sec^ eSs' - p'" " pf^^ tan eBs' sec^ eBsj 



(5.7) 



(5.8) 



(5.9) 



In writing this expression, we have used the notation of p"'! , for example. This signifies a component 
of p which can take only the 'parallel' indices, i.e., and 3, and is moreover different from the index a 
appearing elsewhere in the expression. 



Using now the definition of p' from Eq. (5.4), we can write 



where 



^Ap - + , 



i?ip = ^Apa||/3 p"" taneSs + p'"" taneSs' 



(5.10) 



(5.11) 



and 



R 



(2) 
Ap 



tan eBs 
— tan eBs' 



'£Ap03 +eApa||/3||P""p'^" + EApa,, ,3^ (p""/^ +p°«p'^^ tan^ eSs') 

m^EApoa + eApa||;3||P'"iip"'" + sxpc.,^p^ (p'^iip"'^ tan^ eBs'^ 



(5.12) 
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Obviously, R^^^ is gauge invariant, i.e., k^R^xp 



kPR 



(1) 

Ap 



0. To simplify the other term, we first note 



that the combinations in which the parallel components of p and p' appear in Eq. ( ^.12 ) can be simplified 
by using the following identity: 



£ Apa II /3|| = -eAp03 O, ■ b\\ , 



(5.f3) 



which holds for any two vectors a and h. For the terms involving the transverse compoents, we make an 
important observation. We will be performing the calculations in the rest frame of the medium where 
p ■ u = pq. Thus, the distribution function does not depend on the spatial components of p. In the 
last term of each square bracket of Eq. (5.12), the integral over the transverse components of p has the 
following generic structure: 



(5.14) 



Notice now that 
d 



dpp^ 



g*(p,s)g$(p',s') 



4;f taneSs A +taneSs>"^^')e*(P'^)e*(p''^'). (5.15) 

However, this expression, being a total derivative, should integrate to zero. Thus we obtain that 

taneSs/^ = ~ tan eSs' p'^^ , (5.16) 

where the sign '=' means that the expressions on both sides of it, though not necessarily equal alge- 
braically, yield the same integral. This gives 



tan eBs 



p 



tan eBs + tan eBs 
tan eBs 



tan eBs + tan eBs' 



kP- 



(5.17) 



Using these identities, we can rewrite Eq. (5.12) in the following form: 

tan eBs tan eBs 



R^x} = £Ap03 - p\ ) tan eBs ~ [m^ - p'^) tan eBs' 

'Ap 



£ApQ|j/3i 



where 



R 



(2a) 
Ap 



taneSs taneSs' , , ~ , ^, 
-^^^""^ taneg(. + .0 ^^^^^ "'^ 



^(2h\ ,7 9n ^ ,7 /2n ^ I tane^Bs taneSs' , , 
R^*"^ = (m^ -p|)tane6s - (m^ -p|p)taneSs' ^ , „a (p + P ) ' k\\ 



(5.18) 

(5.19) 
(5.20) 



taneS(s + s') 

The term called i?*^^**) does not vanish on contraction with arbitrary k^ . This term is not gauge invariant, 
and therefore must vanish on integration. In the Appendix, we show that this is indeed true, so that the 
contribution to the vacuum polarization tensor which is odd in B is given by 



(27r)4 



R 



(1) 

Ap 



R 



(2a) 
Ap 



= Me?£\pa.i)k^ 



d-^p 
(27r)4 



77- (p) / dse*(P'^) / ds' e*(f''^') 



^ ^ ^ / taneSs tane;Bs' , , ~„ 
tane6s + p'"« taneSs' ^ (p + p') " 



(5.21) 



tan e;B(s + s') 

In order to perform this integral, we need to introduce further assumptions, which will be done in Sec. ^. 
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6 Dispersion relations 



6.1 Magnetic field-independent terms in the vacuum polarization 

The contributions to the vacuum polarization tensor determines the equation of motion of a photon 
through the medium. We have aheady found the magnetic field-dependent terms in the vacuum polar- 
ization. To obtain the dispersion relations, however, we need also the terms which which do not depend 
on the background magnetic field. These terms are necessarily even in B and therefore did not appear in 
0\p. Here we outline the calculation of these terms. 

Rather than going back to Eq. ([5.6[) which contains also the even terms in B, we use directly the 
propagators at B = given in Eqs. ( |3.15[ ) and ( 3.16| ) to write the background dependent dispersive terms 



as 

^4 



n'Ap(fc) - -^e' J ^ tr [7A5o^(p)7p^5<r (P') + 7a ^S^{ph, S^^ip')] , (6-1) 
Changing, as before, the integration variable in the second term, we obtain 

W (k) - 2 f d^P Sjp^-m^) 
' J (2^)3 {p + ky-m^ 

xtr [7a -I- ■m)r]F{p)lp{i> + ^ + m) + 7p(|5 - m)r]Fi-ph\{^ + ^ - m)] 

(2^ k^ + 2p-k ^"^P^PP +P\kp + kxPp - gxpP ■ k] [/+ + /-] . (6.2) 

In writing the last form, we have put — m? in the denominator and in the trace, in view of the presence 
of the (5-function, and used 

VF{p) + ilF{-p) = f+ + .f-, (6.3) 

where we introduce the notations 

/±-/F(bo|,TAi). (6.4) 



The expression presented in Eq. (3.2) has a particularly simple form in the long wavelength limit, 
i.e., in the limit oi K — Q. In this case, one can show that the IIoo component vanishes, whereas the 
Ily components are proportional to the unit matrix. Since the same is true for the tensor y^j,, we can 
summarize all this information by writing 

nAp(fc) = ^0 5Ap, (6.5) 
where is called the plasma frequency, and is given by 

-0^ - / fl-^) + (6-6) 



{2iTf2Ep \ 3E. 



where P — \p\. 



6.2 Dispersion relations and Faraday rotation 

We have thus obtained expressions for the vacuum polarization tensor. For the rest of this paper, we 



will consider only photon propagation along the direction of the magnetic field. Thus, in Eq. (5.21), 
the index /3 can only take the values or 3. Since the index a appearing in that equation had also 
parallel components only, the antisymmetric tensor now implies that 11"^ vanishes unless both A and p 
are transverse, i.e., have values 1 or 2. Thus, the only non-vanishing components of ^'\p are: 

Il'l^{k) = -U'^,{k) = ~ta, (6.7) 
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where a has to be determined by evaluating the integral in Eq. (5.21). The contributions which come 
from the medium even without the magnetic field has been given in Eq. ( |6.5D . 

To obtain the dispersion relations, we go back to the Lagrangian given in Eq. ( |2.5| ). The equation of 
motion obtained from this Lagrangian is 



In view of the Lorentz gauge condition kpA^ — 0, this can also be written as 



(6.8) 



(6.9) 



For the transverse components of the photon field A^, the above equation implies the following 
condition: 



-la 



A2 



0. 



The eigenvalues of the matrix give the dispersion relations 

fc^ = ± a 

for the normalized eigenmodes 

{A^±iA2)/V2, 

which describe circularly polarized states of the photon. 

Writing fc^ as liP — ^ we obtain the following solutions for K: 

„ li/2 



K± = 



It 



For small magnetic fields, a will be small, and then we can write 

a 



K± 



-ujI 



It 



2(^2 _ ^2) 



which gives, for the difference of the two solutions, 

AK 



(6.10) 

(6.11) 
(6.12) 

(6.13) 

(6.14) 
(6.15) 



For a plane polarized electromagnetic wave propagating with a frequency uj, this means that, after 
travelling a distance £, the plane of propagation will be rotated by an amount t AK. Thus, the rate of 
rotation of the polarization angle is given by 



— = AK = 
M 



(6.16) 



This is the Faraday rotation per unit length. The magnitude of this quantity is thus determined once we 
determine a and wq-Q 

In what follows, we find out what ujq and a are for different types of backgrounds, and consequently 
what is the amount of Faraday rotation sufferred by plane polarized light in such backgrounds. We will 
do this for three different kinds of backgrounds, depending on the relative importance of the temperature 
T — 1/(3, the chemical potential fi, and the electron mass rrie. 

^Instead of a;o, one can also u se the index of refraction r, defined by the relation r = K/w. In absence of the magnetic 
field, i.e., when a = 0, Eq. gives 

We can use this relation to eliminate uio from the formulas above and express everything in terms of r, i.e., the refractive 
index in absence of the magnetic field. 
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7 Results for different backgrounds 

7.1 General observations and assumptions 

Before starting with any of the specific cases, let us note some general features and some common 
assumptions in the calculations. We will perform all the calculations assuming that the background 
medium is at rest in the frame in which we have a purely magnetic field. In other words, for the 4-vector 
u, the only non-zero component is the time component, which has the value unity. All other components 
are zero. 

As already mentioned, we will consider photon propagation along the z-direction (positive or negative). 



In addition, we will take the long wavelength limit, i.e., K <^ co. This implies that in Eq. (5.21), the term 
with the external factor of = uj dominates over the one with = K. 

Finally, we will assume the magnetic field to be small, so that we can use only the linear terms in B. 



To this order, then, the dominant contribution of Eq. ( 5.21 ) is given by 

n'lpik) = 8ie^Bexp3o^I , (7.1) 

where 



s + s' 



Here, since the other factor is already linear in i3, we have used K = in the exponents. Moreover, we 
have made the further assumption that oj <C TOe, which enables us to neglect fco compared to po in the 
factor inside the square bracket. In the notation introduced in Eq. ( |6.7D , we can write 

a = Se^BujI. (7.3) 

The expression for / can be put in a convenenient form. For this, we first define the integral 

^^_{p)pj e-('''-'"')~^l^l / ds' e-'^P"""')-^!^'! r . (7.4) 

(^l") J-oo Jo 



s' - = (s + s') -s' + , (7.5) 



If we now rewrite the factor in the square brackets in Eq. (7.2) in the following form: 

s + s' 
it is easily seen that 

d f 

I ^i-g^Jo- Ji+i J d{m^) J2. (7.6) 
The task is now to evaluate J„ for n = 0, 1, 2. The s-integral in J„ gives 

ds e''*(p'-'"')-^l''l = 27r5(p2 _ m^) , (7.7) 



whereas the s'-integral gives 



Writing now 



[p''^ — m^j 



5 _ ^2) ^ _L [5{p^ _ Ep) + S{po + Ep)] , (7.9) 
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and using 

V-{p) = sign(po)(/+(p) - f-{p)) 
which foUows from the definitions in Eqs. ( 3.12| ) and ( ^.8| ), we obtain 
d*p posign(po) 



J n 



{2n) 



2E„ 



[S{po - E.p) + 5{po + Ej,)] 



U-f- 



{k^ + 2pou;-2PKcos9') 



,\n+l 



(7.10) 



(7.11) 



Here, = (uj,K), P = \p\ and 9' is the angle between K and p. We have denoted this angle by 9' 
in order to emphasize that, for a general direction of propagation, it can be different from the angle 9 
which is measured from the z-axis, i.e., from the direction of the magnetic field which we have already 
specified. For our specific case of photon propagation along the magnetic field direction, however, we will 

put 9' = e. 

We further notice that we can neglect the term fc^ because of our assumptions stated earlier. Thus, 



Jn 



(fp 

(2^ 



(/+ - /-) 



1 



1 



{Epuj - PK cos 9) 



n+l 



{-EpUj - PKcost 



\n+l 



(7.12) 



The azimuthal integration gives a factor 2tt, and the 9 integration can be exactly performed here. This 
shows that J„ = for even values of n. This conclusion can be avoided only if uj and/or K becomes 
comparable to me- Since we have already assumed otherwise, we get 



-Ji, 



(7.13) 



I = - 
8 



d^p 



if+-.f-) 



{Epuo - PK cos 9Y {-EpUj - PK cos 9y 



(7.14) 



In general, however, even this integral cannot be performed analytically. So, in order to discuss the 
amount of Faraday rotation caused by this term, we need to take recourse to some specific limits. 



7.2 Connection with the form factors 

It is of interest to see how our final result for H'^^ conforms to the general form obtained on the basis of 
gauge and Lorentz invariance, a subject that was discussed in Sec. ^. At the end of Sec. ^, we remarked 
that in our case, we can get at most four independent form factors, viz, those associated with the field- 



dependent tensors P'xp\ 



p: 



,//(3) pit{7) 



and P 



"(9) 



However, the simplifying assumptions made above 

"(1) 



imply that all components of u\ vanishes to the leading order. Thus, only P^^ survives in this case. 
Moreover, since we choose the direction of propagation to be along the magnetic field, (k-B)\ = as well. 
Thus, from Eq. (2.17), we find that the tensor P'^p^\ in the case of our choice, is simply proportional to 
B\p. This is what the explicit calculation of Eq. (7.1) tells us as well. 



7.3 A non-relativistic background 

Suppose we have a gas of electrons and positrons where all the particles are non-relativistic. In this case, 
we can put Ep « me within the integral, and neglect all occurrences of P since it is small compared to 
Ep. Then we obtain 



1 



Am-^LO^ 
1 



d^p 



(27r)3 

{Ue - Tie) ■ 



{f+-f-) 



(7.15) 
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Using Eqs. (7.3) and (6.16) now, the Faraday rotation per unit length is obtained to be 

(rie - ng) , 



(7.16) 



where woj in this limit, can be simplified by using the general formula in Eq. (6.6): 



2e2 



(2^ 



[/+ + /-] 



(ne + ng) 



(7.17) 



If the background contains no positrons, the expression for Faraday rotation can be written as 



where lOc = eB/nie is called the cyclotron frequency. 



(7.18) 



7.4 A degenerate background 

We now consider a degenerate electron background at zero temperature. The distribution functions are 
now given by 



/+ 
f- 



for P <Pf 
for P> Pf 



0. 



(7.19) 



where Pp is called the Fermi momentum. As we know, although the temperature is zero, the electrons 
need not be non-relativistic in this case, since Pauli exclusion principle would require all electrons to be 
in different states, and so some of them can be at very large momentum. The number density of electrons 
in this case is given by 



rifi = 2 



(2^ 



n3 

37r2 



(7.20) 



In this case, we first calculate the plasma frequency. Performing the angular integrations of Eq. ( |6.6D , 
this can be written in the form 



dx 



(l+a;2)l/2 3(1+2.2)3/2 



(7.21) 



where x is the integration variable defined by P/nie, and xp = Pf I'nie- The integration can be performed 
in a straight forward manner by substituting x = sinh C, and the result is 



e^mg x\ e? Pp 

yrT4 " 3^ ^ 



where Ep \s the Fermi energy. 



Ep = X P\ 



(7.22) 



(7.23) 



We now evaluate the integral /. Starting from the expression in Eq. (7.14) for the general case, we 
perform the angular integrations to obtain 



/ = 



16 

8^2, 



1 r°° 

-2-/ dpp{,u-f^] 
1 r°° p2 



1 



1 



Epuj - PK Epuj + PK 



(7.24) 
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Using now the distribution functions appropriate for this case from Eq. ( |7.19D , we obtain 

I = o I 9 Pf - TOetan~^(PF/TOe) , 

where the result of the arctan function is restricted within the domain to tt/2. From this, we obtain 
the Faraday rotation per unit length to be 



p3 



Pf 



TOg tan ^{Pp/mf,) 



(7.25) 



(7.26) 



It can be easily checked that if Pp <C me, in which case the background is non-relativistic, the formulas 



derived for this case reduce to those derived in Sec. 7.2 



7.5 An ultra-relativistic background 

Let us now discuss the case where the temperature T of the background is much higher than the electron 
mass. In this case, we can put Ep « P. Then, using the dimensionless integration variable y = P/T, the 
plasma frequency can be expressed as 

"'0 ^ 3^ i y Uxp(y~/3A.) + 1 + exp(y + /3;.) + lj ' ^^'^^^ 

This integral can in fact be performed exactly. In the first integral, use the new integration variable 
y' — y — Pfi. In the second one, use y' = y + PfJ.. The resulting integrations can then be written in the 
form 



2e2 



dy' 



22/' 



ef' + 1 



dy 



/3m 



ey' + 1 



dy' , , 



(7.28) 



The first integral can now be performed by expressing the denominator as a geometric series. The other 
two can be combined after substituting y' — — y' in the second integral, and the final result is 



1 



2 1 



37r2 



(7.29) 



For the integral /, we start from the general expression in Eq. (7.24). Using similar substitutions as 
before, we obtain 



/ = 



1 



2,, ,2 



dy 



1 



1 



exp(y - /3/i) + 1 exp(y + /3/i) + 1 



87r2w 



The Faraday rotation is obtained from Eq. ( 3.1(: ) and (7.3): 



(7.30) 



(7.31) 



If we put (3 — > oo, the results obtained for this case reduce to those obtained in § 7.4 with TOg = 



As for the previous cases, one may want to express these results in terms of the number densities of 
electrons and positrons in the medium rather than in terms of the chemical potential fj,. The connection 
comes from the relation 



rip 



(27r)3 



(/+ - /-) 



(7.32) 



Again, the integration can be performed exactly, following the steps described above, and the result is 

Tip — = + ^r—^ . (7.33) 



3/32 



ill 
37r2 



One can use this to express /i in terms of rig 



ng. 
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8 Conclusions 



We have thus shown that the amount of Faraday rotation depends very significantly on the characteristics 
of the medium in which the background magnetic field rests. If the medium consists of non-relativistic 
particles only, we obtain the formula given in Eq. ( 7.18 ). Usually, we assume that this formula should be 
applicable for low-temperature media, since the particles should be non-relativistic in this case. However, 
we show that if the medium is strongly degenerate, the formula changes. For Pp <C me it still agrees 
with the non-relativistic result as it should. But for Pp 3> me the change is drastic, and Faraday rotation 
becomes very small, as can be seen from Eq. (7.26). Similarly, if the medium is so hot that the kinetic 
energies of the particles are much larger than their masses, we obtain a different result, as shown in Eq. 
(7.31). But it is interesting to note that in all the cases discussed, the quantity a has the same dependence 
on uj, viz., that it is inversely proportional to uj^/uP^—lJ^. 



Appendix: Proof of gauge invariance 



In the text, we claimed that the contribution coming from R^"^^^ must vanish in order that the vacuum 
polarization tensor is gauge invariant. Here, we justify this claim. This contribution is proportional to 
the following integral: 



C = 



(27r 



^T^r J-oc Jo 

I 1 2\ . I 1 i2\ . rj I tau eSs tau gSs' , ,, , 

(m^ - pi)ta.neBs - (m^ ~ p.ntaneBs — — (p + p ) ■ kn 

" " taneo(s-fs') " 



Using the definition of the exponential factor <i>(p, s) from Eq. (3.6), we notice that 

taneSs e*'^'^' e*(p''^') = taneBs [i-^ + (pf - sec^ eBs' p'l)} e*'"'^) e*(f''^') 



m2tanees'e*(P'^)e*(p''"') = tan eBs' {i— + {pj - sec^ eBs pi)} e 
This implies that we can write 



C = Ci+iC2., 



where 



j77_(p) / dse^^P^''> / ds' e^^P'-^'^ 



d^p 

W) 



{pf - sec^ eBs' p'1 - p|) tan eBs - (p^ - sec^ eBs p\ - p'^) tan eBs' 



taneSs taneSs' , 
taneS(. + .') ^^^-^^h 



Co 



f d'^p |-oo |.oo / ^ 

/ . N 1 7?- (p) / ds ds' { taneBs -;— ~ ta,neBi 
J (27r)4 J_^ 7o V ds' 



:'± 
ds 



(A.l) 

(A.2) 
(A.3) 

(A.4) 



(A.5) 



(A.6) 



Let us first consider the contribution C2. Performing the s'-integration in the first term and the 
s-integration in the second, it can be written as 



C2 



d-^p 



V- (p) 



a'^ip',s') 



ds tane6se*(P'") _ e*(P'^) 

) 
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ds' taneBs' e*(P''"') 



(A.7) 



The second term vanishes since e**^^'*^ vanishes at both hmits due to the term — e|s| in it. The other 
exponential survives only at the limit s' ~ 0, and gives 

C2 



^rj^{p)J tane6se*(f'^) 



- 0, 



(A. 



where the last step follows on performing the integration over p, since ^(p, s) is an even function of p 
and 77_ (p) is odd. 

Let us now look at the other contribution, Ci. Separating out the terms involving parallel components 
from those involving transverse components, we write 

/j4 poo poo 

J^V-{P) J ds ds' exp[$(p,s) + $(p',s')] 



(pf ^ P\\)[ '^^^ + ^ 



tan eBs tan eBs' 



taneS(s + s') 
+p\ tan eBs' sec^ eBs — p']_ tan eBs sec^ eBs' . 

From the definition of $, it follows that, apart from the small convergence factors, 
<^{p,s)^<^{p',s') = l(s + .')(pf +P|-2m2)-l(.-.')(pf 



(A.9) 



~7b 

i 

eB 



'tan eBs' p'l + tan eBs p\) 

(p'^ +p\- 2m2) ^ _ (pj2 _ pjj) ( _ p'l tan(C ~ ~ pI tan(C + C) 



(A.IO) 



eB[s + s'). 
eB{s-s') 



where we have defined two new parameters 

f = 

c = 

Thus, 

ieB Ae*(p.-)+*(p'.-') = e*(p,-)+*(p',-') (pj2 _^2_ ^^^2^^ _ ^ ^^^2^^ _^ 
Using this, we can rewrite Eq. ( [A.9| ) as 

X ftane&s + taneS.'- taneg. tanegsV^^ 

V tanee(s + s') / dC 

+e*(P,s)+*(p',.') I tan eBs' sec^ eSs' f 1 ^^i^^ 

1^-^ V taneS(s + s') 

2 ^ K> 2 » /"i taneSs' 

— p I taneossec eos 1 — — 7 

V tan eS(s + s') 



(A.ll) 



(A.12) 



(A.13) 
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We are now left with only the transverse components everywhere except in the exponents. To write 
them in a useful form, we turn to Eq. (5.15) and take another derivative with respect to p"^. From the 
fact that this derivative should also vanish on p integration, we find 



a 8 ° 
P±P'1 = 



1 



tan eBs + tan eBs' 



ieB 



tan^ eBs' 



tan eBs + tan eBs 



In particular, then. 



1 



tan eBs + tan eBs' 



It then simply follows that 



p1^ 



1 



tan eBs + tan eBs' 



ieB- 



-ieB 



tan-^ eBs' 
tan eBs + tan eBs' 



tan^ eBs 
tan eBs + tan eBs' 



(A.14) 



(A.15) 



(A.16) 



We now put these into Eq. ( A.13| ). After some straight forward but cumbersome algebra, it is found that 
the terms involving k cancel out, and we are left with 



Ci = ieB 



iv-ip) 



(2^)4 
tan eBs + tan eBs' 

g*(p,s)+<I'(p',s') 



ds 



ds' 



t'ATieBs tan eBs' \ d 
taneS(s + s') ) dC, 



^<S>(p,s)+<t>{p' ,s') 



tan eBs + tan eBs' 



tan eBs' sec^ eBs' ( 1 - 



tan eBs sec^ eBs {l 



tan eBs 
tan eB{s + s'] 

tan eBs' 



taneS(s + s') 

It is straight forward to show that this can be written in the following form: 

d^p 



ieB 



(27r) 



/ ds ds'-T{tO. 
J^oo Jo dC 



where 



— tan eBs + tan eBs' — 



tan eBs tan eBs' 
taneS(s + s') 



^*(p,s)+*(p',s') 



(A.17) 



(A.18) 



(A.19) 



with s and s' related to ^ and ^ through Eq. (A. 11 



We can now change the integration variables to f and This gives 



d^p 



eB J (2^)4 I / ® " -H^^ 



d 

dC' 



dC 



d 
dC 



2i 



d*p 



eB J {2ttY 



(A.20) 



since the other term vanishes at the limits. In this integrand, C — which means s' = 0. Looking back 
at the definition of J-, we find 



eB' 



tan 2^ . 



(A.21) 



This is an even function of p, whereas {jp) is odd. Thus, the expression vanishes on integrating over p. 
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